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■ Abstract 

We explore the usefulness of the existing relations between the 
S'-matrix and time delay in characterizing baryon resonances in pion- 
. nucleon scattering. We draw attention to the fact that the existence 

J> I of a positive maximum in time delay is a necessary criterion for the 

existence of a resonance and should be used as a constraint in con- 
ventional analyses which locate resonances from poles of the S'-matrix 
00 , and Argand diagrams. The usefulness of the time delay plots of reso- 

nances is demonstrated through a detailed analysis of the time delay 
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■ in several partial waves of ttN elastic scattering. 
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• 1 Introduction 

: 

^ I Ever since the discovery of the first excited nucleon state [1], the baryon 

^ I resonances have played a major role in particle and nuclear physics and have 

contributed crucially to the search of the fundamental building blocks of 
nature. We perceive them now as the low energy manifestation of quan- 
tum chromodynamics (QCD) with three quark degrees of freedom. However, 
low energy QCD is still not well understood and very often one is left with 
models. The status of the resonances can differ from case to case as can 
their parameters extracted from different experiments. The N* programs at 
Jefferson Lab [2] and the forthcoming Japan Hadron Facility [3] have both 
revived the area. The various theoretical studies [4-7] and the hope to find 
exotic hadrons [8], make the area once again an exciting field of physics. It 
is then justified to look at the baryon resonances from a yet different per- 
spective and to analyze the existing data in a novel albeit well established 
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way, while waiting for new experimental results. Specifically, we refer to 
the time delay method which was introduced into scattering theory and es- 
pecially resonance physics by Eisenbud and Wigner [9, 10, 11]. In view of 
the considerable number of standard textbooks [12-17] and numerous papers 
[9-11,18-39] published since the seminal paper by Wigner [10], we provide 
only a short introduction. Time delay is a measure of the collision time in 
a scattering reaction which can be calculated directly from the phase shift 
or the T matrix. Obviously, such a concept has a close connection to the 
appearance of an unstable intermediate state (resonance), which, due to its 
finite lifetime, "delays" the reaction. Though the interest in time delay ever 
since the first papers was unabated, it is only recently that it has been used 
in practice in quantum scattering theory (in chaotic scattering [40], hadron 
resonances [33,37-39] and heavy ion collisions [35]) and tunneling phenom- 
ena [36], with success. The present work is simply a logical extension of this 
program carried over to baryon resonances, partly done already in [37]. 

To identify the baryon resonances, one performs a partial wave analy- 
sis of the meson baryon scattering data and obtains the energy dependent 
amplitude (or T-matrix) by fitting cross section data. Resonances are then 
determined by locating the poles of the T-matrix on the unphysical sheet 
and studying the Argand diagrams of the complex T-matrix. Due to model 
dependence in the analyses of the energy-dependent amplitudes, there are 
differences in the resonance parameters quoted by different groups [41]. The 
resonance receiving confirmation from several analyses is considered to be 
well estabhshed. Though we do not dispute the usefulness of the pole of the 
S'-matrix, we note that there exist several views in literature, regarding the 
definition of a resonance. In a review article [42], Dalitz discussed various 
criteria for the existence of a resonance elaborately, with the conclusion that 
for the case of a pole in the S'-matrix, S{E), in the unphysical E-plane lying 
sufficiently close to the physical E axis, there is no ambiguity in the conclu- 
sion of the existence of a resonance. However, the authors in [43] constructed 
examples in such a way that a sharp resonance was produced without an ac- 
companying pole in the unphysical sheet. They noted that even the inverse 
correspondence, namely, (pole of the S'-matrix on the unphysical sheet) 
(unstable particle) may be questioned. In [44], it was pointed out that a 
peak in the cross section cannot be conclusive evidence of a resonance. In 
[45], in addition to time delay, the exponential decay law was required as a 
signal of a genuine resonance (this may be in view of the existence of double 
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poles, which would lead to a non-exponential decay [46]). Cautious remarks 
on the use of Argand diagrams can be found in [19, 47]. The many different 
opinions reflect only the fact that the issue is not yet satisfactorily settled. 
Indeed, unstable particles remained to be problematic even until now [48]. 
We make use of the requirement stated in literature and text books [9-39], 
namely, the formation of a resonance should introduce a large positive time 
delay in the scattering of particles. We try to extract resonance parameters 
from the energy distribution of time delay by locating the position of the 
local maximum and reading off the width as advocated e.g. in [49]. Though 
the non-resonant background can deform the positive resonant structure in 
the vicinity of a resonance, we do expect some positive region around the 
resonance point, with perhaps a less dominant peak. This is confirmed by 
our study. 

Starting with the definition of time delay in terms of the S'-matrix, we 
obtain its relation with the T-matrix and scattering phase shifts. We shall 
first demonstrate the usefulness of the method with examples of well-known 
N and A resonances. Later on we proceed to the analysis of time delay 
in various partial waves of ttN elastic scattering, using the available single 
energy values as well as some energy-dependent forms of the T-matrix. Before 
we move on to the discussion of time delay, it is important to note that the 
time delay plots of the present work are not the same as speed plots [50] which 
have been sometimes referred to as time delay plots in literature. Speed plots 
are positive definite by definition. Time delay plots can also assume negative 
values and only a positive peak signals a resonance. In the elastic region, the 
speed is equal to time delay up to a constant factor, but once the inelastic 
channels open up, this is no longer true [37]. 

Considering the fact that the time delay method has so far not been 
applied to baryon resonances (but has been successfully applied to meson 
resonances [39]), our study is a practical test of time delay when apphed 
directly to data. When applied to theoretical T matrix solutions, we could 
say that indeed the model is being tested, if we consider the resonances to 
be well established. 

In passing, we note that time delay is also related to the so called arrival 
time in quantum mechanics [51] and has also been used to obtain the density 
of resonances [52] . 
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2 Time delay in resonant scattering 



We shall now discuss the expressions which quantify time delay and can hence 
be used to characterize resonances. 

2.1 Relation to phase shifts 

In the early fifties, using a wave packet analysis, Bohm [15], Eisenbud [9] 
and Wigner [10], obtained an expression for the time delay At in binary 
collisions. In the case of elastic scattering, they derived At in terms of the 
energy derivative of the scattering phase shift as follows: 

At = 2S-, (1) 

The formation of a resonance in a scattering process, introduces a positive 
time delay between the arrival of the incident wave packet and its departure 
from the collision region. From the above relation, one expects the phase 
shift to increase rapidly in the vicinity of a resonance. 

The wave packet analysis of time delay was extended by Eisenbud to 
inelastic collisions [9]. He defined the delay time matrix At, such that an 
element Atij of this matrix, corresponded to the time interval between the 
outgoing wave in channel j and the ingoing wave in channel i. This time 
delay, At^, is related to the ^'-matrix as follows: 

At,,^Re[-ih{S,,r'^]. (2) 

Before we proceed further, we note that the phase shifts, in principle, 
depend on the orbital angular momenta, /, of the initial and final states 
respectively and on the total angular momentum J. However, we have sup- 
pressed this dependence in the expressions whenever not relevant. In the 
present work, we consider ttN elastic scattering, which is the scattering of 
a spin zero and spin one half particle. Since the total spin in the final and 
initial state is S = S' = 1/2 and conservation of parity gives / = the total 
angular momentum J takes the values 1 — 1/2 and 1 + 1/2. The ^'-matrix is 
diagonal in I and its elements are related to phase shifts as Sfi — exp{2iS'^), 
for the elastic case in the absence of inelasticities. 
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We see that in the case of purely elastic scattering (j = i), and using a 
phase shift formulation for the ^'-matrix where S = e^^^, we get, 

Atu = 2h— , (3) 

which is the same as Eq. (1). These Atu are related to the lifetimes of 
metastable states or resonances in elastic scattering (see [18]). At high en- 
ergies, where apart from elastic scattering, the possibility of scattering into 
inelastic channels also opens up, the elastic (S-matrix element is defined as 
S = rye^"^, where r] is the inelasticity parameter defined such that < 77 < 1. 
Substituting the modified S (i.e. S — rje^^^) in Eq. (2), gives. 



At„ = Re 



, ^. dS df] 1 



The above equation is the same as Eqs (1) and (3). Thus it can be seen 
that the expression for the time delay, Atu, for elastic scattering is the same, 
irrespective of the presence of inelastic channels. 

It is clear from the above expressions that time delay can also take nega- 
tive values resulting from phase shifts which decrease as a function of energy. 
However, the negative delay times cannot assume arbitrarily large values. In 
the case of elastic scattering (for the case of Z = and 1) it was shown by 
Wigner [10], that the causality condition puts a constraint on the lower value 
of the phase shift derivative (related in an obvious way to time delay), which 
in case of high momenta, i.e. for large k is given as, dSi/dk > —a. a can 
be interpreted as the range of the interaction potential. We do observe some 
regions of large negative At which will be discussed in Section 3. 



2.2 Relation to T-matrix 

Instead of using the phase shift formulation of the ^'-matrix, we now start 
by defining the S'-matrix in terms of the T-matrix, i.e., 

S = l + 2iT, (5) 

as is usually done in partial wave analyses of resonances [53, 54] . The matrix 
T contains the entire information of the resonant and non-resonant scattering 
and is complex (T = T^+iT^). Substituting from Eq. (5) into the expression 
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for time delay in (2), the time delay Atjj, in terms of the real and imaginary 
parts of the amplitude T is given as, 



+ 27;f— ^ -2T^ 



" " dE 



(6) 



where S'*jS'jj can be evaluated using Eq. (5). In the present work, we have 
evaluated the time delay in TriV elastic scattering and hence, i corresponds 
to TiN in the above equation. 

Although a simple Breit-Wigner (BW) is not always a good choice to 
describe a broad hadronic resonance, it is instructive to see the results we 
get for time delay, starting from a BW matrix element. If we insert one such 
commonly used form of the T-matrix [54] in resonance regions, namely, 

T= p) 

£„-£-ir/2' 

in (6), we obtain, 

and the time delay at the resonance energy Er (within the assumption that 
the widths are not energy dependent) is, 

At{ER)Bw = Y- (9) 

A simple BW T-matrix, however, can be misleading, especially while dis- 
cussing time delay. The reason among others is that it lacks certain usually 
expected properties (threshold behaviour being one of them). We shall come 
to this point in greater detail in section 4. 

Before ending this section, we note the dependence of time delay on wave 
packets. It is well-known that the survival probability and lifetime of an un- 
stable quantum state depend on its preparation. Explicit formulae including 
wave packets can be found for unstable neutral kaons in [55]. Wc expect a 
similar dependence to be present in the expressions for time delay. Indeed, 
as given in [16], 

/■oo d^ 

AtiE) = Sn'h dE'\A{E',E)\'2^ , (10) 
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where A{E') is the initial wave packet in momentum space. If the wave packet 
is sharply centered around an energy E, we recover Eq. (1). In scattering 
processes where one measures the cross sections and distributions, the wave 
packets are indeed narrow, i.e., the energy spread AE <^ F (see the second 
and last reference in [13] for a discussion of this issue). Hence we can use 
Eqs (1-4) to calculate time delay. 

In the next section, we shall evaluate the time delay in several partial 
waves of nN elastic scattering. We have checked that the values of time 
delay, Ata, obtained either using the derivative of the real phase shifts as in 
Eq. (3) or the T-matrix as in Eq. (6) arc the same. Since both the methods 
are equivalent, one can in fact use fits to the single energy values ^ of phase 
shifts to extract resonance parameters. 

3 Time delay plots of resonances in ttN elastic 
scattering 

We now analyze the existing nN scattering data using time delay plots. To 
demonstrate the usefulness of the method, we plot time delay in the energy 
regions where two well-known baryon resonances occur. In Fig. 1 are shown 
the real and imaginary parts of the complex T-matrices, the phase shifts 
and the corresponding time delay in the P33 and D13 partial waves in ttN 
scattering, evaluated using the T-matrices (solid lines) which fit the single 
energy values of T very well. The filled circles in Fig. 1 are the single energy 
values of phase shifts extracted from the cross section data on ttN elastic 
scattering [56]. The widths of the P33 and i^i3 peaks at half maximum can be 
read from Fig. 1 to be around 116 and 50 MeV respectively. The peaks in the 
energy distributions occur at 1216 and 1512 MeV respectively. The average 
values of Breit-Wigner masses (widths) given in the Summary Table (ST) of 
the Particle Data Group [41] for these P33 and D13 resonances arc 1232 (120) 
and 1520 (120) MeV respectively. The A(1232) decays almost 100% to the 

-"^Thc values of phase shifts in different partial waves obtained by fitting the cross 
section data at the available energies are known as single energy (SE) values of phase 
shifts. The error bars on these phase shifts naturally depend on the errors in the measured 
cross sections. The elastic T-matrix element is related to the phase shift and inelasticity 
parameter, rji, as: Ti = {rjie"^^^' — l)/2i. Thus, one can also obtain SE values of the 
T-matrix. 
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Figure 1: Single energy values of the real part of the T- matrix (filled tri- 
angles), imaginary part of T (open squares), phase shifts (filled circles) and 
the time delay At evaluated in the P33 and D13 partial waves of nN elastic 
scattering. The time delay is evaluated using the T-matrix given by the solid 
lines which fit the single energy values very well. 

ttN channel and hence the time delay width seems to be in good agreement 
with the above value listed in the ST. The D^^ has a branching ratio of 50 
to 60% to the ttN channel and the width of the time delay distribution is 
consistent with the partial width listed in the ST. Thus we see that in the 
case of purely elastic scattering as well as in the case of elastic scattering 
in the presence of inelastic channels, the method is quite useful. The peak 
position and width of the time delay distribution give the mass and elastic 
partial width of the resonance, respectively. 

Interestingly, the P33 phase shift of the only nN resonance (A(1232)) 
in the elastic region, remains positive and shows the characteristic resonant 
jump in this region. Hence, in this case, the speed defined in [50] is the same 
as time delay up to a constant factor. 
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3.1 New resonances from single energy values of phase 
shifts 

We shall now evaluate time delay from fits to single energy (SE) values of 
phase shifts. Since the results depend crucially on the quality of the data, we 
chose data sets with small error bars and made separate n^^ order polynomial 
fits to different energy regions of the phase shift. It would be more appro- 
priate to consider error bars and perform a fit, with a certain function. 
However, such a procedure would not be able to pick up the small structures 
and would amount to giving results similar to the energy dependent ones. 
We also chose to fit SE values of phase shifts rather than the SE values of 
real and imaginary parts of the T-matrix, simply as a matter of convenience. 
The time delay evaluated using fits to phase shifts or T-matrices is actually 
the same. The advantage of calculating time delay from such fits is that the 
results are directly related to data. The disadvantage is that they are sensi- 
tive to the quality of the data and hence to the fit. There also exists the well 
known continuum ambiguity problem with the SE values of phase shifts [58]. 
However, the present work does not aim at finding solutions to the problems 
related to the extraction of SE values. Hence, we use the values as available 
in literature and check if we still get some useful results for time delay. 

We perform this analysis for the / = 1/2 partial waves, Pn, P13, D13, 
Sii and Fi5 in t^N elastic scattering. We note that in spite of the above 
mentioned problems, we get strikingly similar peak positions and widths as 
compared to the Summary Table resonance parameters. 

The results in Figs 2 and 3 reveal that time delay has the following main 
characteristics: (i) it locates well established resonances (ii) the positive 
peaks are more prominent than in the case where we calculated the same 
quantity for the energy-dependent solutions (see section 3.2 below), (iii) there 
exist regions of negative time delay in addition to the positive peaks (iv) the 
new feature here is that we find additional resonant peaks. At present, given 
the quality of the data, it is not clear if these new structures are artifacts of 
the fit to the data or genuine indications of (new) resonances. For example, 
in the Pn and Su cases, we have hints for new resonances in the higher 
energy regions where the quality of data is worse. On the other hand, some 
one-star resonances like Pii(2100) and 5'ii(2090) are not excluded. We find 
evidence for the 3-star resonances Pi5(2000) and Pi3(1900), with some indi- 
cation that the latter consists actually of two nearby resonances. In the D12, 
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Figure 2: Single energy values of phase shifts and the corresponding time 
delay in the Su partial wave evaluated from a smooth fit to the phase shifts. 

partial wave, we observe distinct peaks at 1512, 1695 and 1940 McV, which 
could be associated with the 4-, 3- and 2-star resonances A^(1520), A^(1700) 
and A^(2080) respectively. Note however that the existence of the two small 
peaks in this case depends crucially on two data points, and hence on the 
fit. These two points are sufficiently above continuum to justify the peaks 
(more so as they can be associated with known resonances). There seems 
to be more structure in the 1400 — 1700 McV region of 5*11. A fit made to 
this detailed structure reveals the possibility of four resonances around 1650 
MeV. Indeed, there is some support for this structure from recent works in 
literature [5, 59], where the existence of new resonances at 1.6 and 1.7 GeV 
is predicted within quark models. 

With the availability of more precise data on cross sections which would 
enable a better extraction of the SE values of phase shifts, one could locate 
the resonances from time delay plots more accurately. We limit the discussion 
in this section only to the 5 partial waves shown in Figs 2 and 3, since a 
detailed analysis with all partial waves would make sense only when the SE 
values of phase shifts would be better known. We list our findings in Table 
1. 
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Figure 3: Single energy values of phase shifts and the corresponding time 
delay in the Pn, Pi,3, Di^ and F15 partial waves of ttN elastic scattering. 
The time delay is evaluated using the solid lines which are fits to the single 
energy values of phase shifts. 
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Table 1. Peak positions and corresponding widths from time delay plots 
using fits to single energy values of phase shifts. Masses and widths are in 
MeV. 



L27,2J 

Mass"*"*"" (Full 
Width) 


ST average 
Pole (P) 
Re P [-2 Im P] 


B.F.=^,^^/^ 
(ST average) 


Partial width 
= B.F.x 
(-2 Im P) 


Time Delay 
Peak [Partial 
Widths] 


Di3 

1520**** [120] 
1700*** [100] 
2080** [-] 


1510 [115] 
1680 [100] 
1824 - 2120 


0.5 - 0.6 
0.05 - 0.15 


58 - 69 
5 - 15 


1512 [48] 
1695 [12] 
1940 [18] 


Fl5 

1680**** [130] 
2000** [-] 


1670 [120] 


0.6 - 0.7 


72 - 84 


1685 [91] 
1940 [33] 


Pii 

1440**** [350] 
1710*** [100] 

2100* [-] 


1365 [210] 
1720 [230] 


0.6- 0.7 
0.1 - 0.2 


126 - 147 
23 - 46 


1440 [207] 
1700 [37] 
1830 [65] 
1940 [13] 


Pl3 

1900** [-] 








1520 [101] 
1975 [51] 
2140 [54] 


Sii 

1535**** [150] 
1650**** [150] 

2090* [-] 


1505 [170] 
1660 [160] 

1795 - 2220 


0.35 - 0.55 
0.6 - 0.8 


60 - 94 
96 - 128 


1510 [37.8] 
1590 [25.2] 

1630 [~ 40] 
1680 [~ 24] 
1700 [25.6] 
1810 [31] 
1920 [38] 
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3.2 Time delay from energy- dependent amplitudes 

One of the standard methods of characterizing resonances involves locating 
the poles of an energy dependent T-matrix on the unphysical sheet. If these 
poles correspond to resonances, a positive maximum in time delay at the 
energies where the poles occur is expected. However, none of the existing 
analyses are constrained by this necessary condition. In what follows, we 
use the energy-dependent solutions obtained from the SAID program [56] 
as an example of such analyses, to evaluate time delay. We start with the 
1 — 1/2 partial waves in nN elastic scattering. In Fig. 4a, we plot the SAID 
solution FA02 of the complex amplitudes. The corresponding time delay, 
using these T-matrices and those from an earlier analysis (SM95 solutions) 
by the same group, is plotted in Fig. 4b. The two solutions give rise to similar 
values of time delay in all except the Pu and P13 partial waves, where the 
peak positions differ. The FA02 solutions which are in better agreement 
with the SE values as compared to SM95 were obtained [60] using a much 
bigger database. The Pii( 2000) peak obtained using SM95 is not seen with 
FA02. Indeed it was noted in [60], that the most significant shifts in the 
pole values occur in the P-waves (Pu and P13). We observe the resonant 
peaks in the Pu, ^n, D12, -D15, P15 and P13 partial waves close to the pole 
positions predicted by the T-matrices. However, small peaks in the Gn and 
Hig partial waves appear at much lower values than the poles. The shifts 
in the time delay peaks as compared to the pole values could be due to the 
presence of a non-resonant background in these partial waves. Yet another 
explanation is offered in section 4. We give a more detailed discussion of the 
resonances in various partial waves now. 

Pii-resonances: The SM95 solution gives a broad peak around 1400—1500 
MeV and a prominent peak at ~ 2050 MeV which could be attributed to the 
Pii(1440) and the less estabhshed Pii(2100) respectively. With the FA02, 
the two peaks at 1370 MeV and 1745 MeV are a clear signal of Pii(1440) 
and Pii(1710) listed in the ST. The FA02 T-matrix has two closely located 
poles at 1357 and 1386 MeV. The broad time delay peak around 1370 could 
actually be due to two closely overlapping resonances. The observation of 
the time delay peak at 1370 MeV which is much lower than the ST value of 
1440 MeV is similar to the finding of Ref. [57]. In fact, most of the time 
delay peaks being close to the pole positions, occur at lower values than the 
parameters of the ST. The prominent peak at 2050 MeV with SM95 is not 
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present in FA02 anymore. 

We also note that Pn resonances at 1500 and 1700 MeV were found in 
[61] from fits to the energy dependence of the amphtude obtained from an 
older VPI single energy analysis. 

Pi3-resonances: We see peaks at 1585 and 1600 McV corresponding to 
the FA02 and SM95 solutions respectively. There exists a pole of the SM95 
solution at 1717 MeV which is close to the resonance Pi3(1720). However 
the FA02 pole occurs at a much lower value of 1584 MeV. Though the pole 
values of the two solutions differ a lot, the time delay peaks with the two 
solutions are quite close. 

Sii-resonances: We observe a positive peak around 1494 MeV, which can 
be attributed to 5*11(1535), again as in the case of Pii(1440) at a consider- 
ably lower value. The positive peak at 1650 MeV is a nice manifestation of 
5*11(1650). We note another phenomenon which commonly occurs in time 
delay plots now. The small peak at 1494 MeV is followed by a large negative 
region (which, as emphasized before, cannot be associated with resonance 
formation). The negative time delay is associated with the opening up of 
new channels in ttN scattering and in fact, the minimum in the dip occurs at 
the energy corresponding to maximum inelasticity. The energy derivative of 
the phase shift (and hence time delay) is related through the Beth-Uhlcnbcck 
formula to the change in density of states in the presence of interaction [52]. 
The negative time delay corresponds to situations where due to the inter- 
action, the density of states is less than in the absence of interaction. This 
is exactly what happens when the inelastic channels open up and there is a 
loss of flux from the elastic channel, thus reducing the density of states due 
to interaction. A detailed discussion on this issue can be found in [37]. A 
repulsive non-resonant background could make an additional contribution to 
the negative time delay. 

Di3 and Fis-resonances: Time delay reveals the £)i3(1520) and Fi5(1680) 
resonances. The three- and two-star resonances Dis^lYOO) and Fi5(2000) 
respectively, do not appear in the time delay evaluated with the SM95 or 
FA02 solutions. This is consistent with the fact that the same solutions do 
not locate the corresponding poles. However, these two resonances do appear 
in the time delay plots obtained from fits to the single energy values of the 
T-matrix. 

Di5 resonance: Though the peak is not prominent, its position and width 
are close to the pole values. 
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Figure 4: (a) Real (dashed lines) and imaginary (solid lines) parts of the 
T-matrix solutions FA02 [56] for isospin, / = 1/2, partial waves, (b) The 
time delay (sohd hnes) evaluated using the T-matrix solutions FA02 (shown 
in (a)) and an earher version SM95 (dashed hnes) by the same group. 
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Table 2. Comparison of nucleoli resonance parameters from time delay eval- 
uated using the FA02 T-matrix solution, with the pole positions [62] of the 
same T-matrix, Summary Table values and Speed Plot poles (P = E - iV/2) 



L27,2J 

Mass (Full - 
Width) 


Speed Plot 
Pole(P) 

T) T) r O T "Dl 

Ke P [-2 Im PJ 


FA02 
Pole (P) 
Ke P [-2 Im P] 


Branching 

to ttN 
decay mode 


Time delay 
Peak [Partial 
- Width] 


^11 
144U (350) 


1385 [164] 


1357 [162] 
1386 [170] 


60 - 70 % 


1370 [298] 
t 


^13 

1520 (120) 


"1 "1 ["1 or\l 

1510 [120] 


"1 "1 /I Ti r\ol 

1514 [103] 


50 - 60 % 


1510 [50] 


■^11 
1535 (150) 


1487 [ - ] 


1516 [123] 


35 - 55 70 


-i A r\ A [ A 

1494 [48] 


■^11 

1650 (150) 


1670 [163] 


1639 [155] 


55 - 90 70 


1650 [145] 


Dl5 

1675 (150) 


1656 [126] 


"1/^/^/1 F"! /I'll 

1664 [141] 


40 - 50 % 


i "1 rr\ol 

1610 [93] 


^15 

1680 (130) 


lb/ 3 [i3oJ 


iO M [i2ij 

1778 [215] 


bU - lO 70 


Ib/U [/oj 


^13 

1700 (100) 


1700 [120] 




5 - 15 % 




^11 
1710 (100) 


1690 [200] 




10 - 20 % 


1745 [50] 


■■^13 

1720 (150) 


1686 [187] 


1584 [287] 


10 - 20 % 


1585 [124] 


X3** 

2100 ( - ) 




2009 [458] 






**** 

2190 (450) 


2042 [482] 


2084 [453] 


10 - 20 % 


1900 [310] 


TJ**** 

^19 

2220 (400) 


2135 [400] 


2230 [553] 


10 - 20 % 


2050 [276] 



t Width is large due to the possibility of 2 overlapping resonances. 
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Gi7 and Hig-resonanccs: The peak positions in time delay are at much 
lower values as compared to the poles and much broader than the partial 
widths corresponding to the poles. It is not clear if such a large shift in these 
2 cases should be attributed to the non-resonant background, since most of 
the other resonances were shifted from the pole values by only few tens of 
McV at the most. Although a remote possibility, it could be that the shifted 
peaks represent resonances without corresponding poles and the poles at 
higher values do not correspond to resonances. These cases could be realistic 
examples similar to the constructed ones in [43] , where the authors doubt the 
one-to-one correspondence between an unstable state and S'-matrix pole. A 
good knowledge of the non-resonant background would clarify the situation. 
An alternative explanation will be given in section 4. 

A comparison of the time delay peaks (evaluated using the FA02 solution) 
with the ST vahics and pole positions of FA02 and Speed Plots is given in 
Table 2. The branching fractions of the various resonances to the vrA^ decay 
channel are also listed. The widths in time delay arc the partial widths 
corresponding to the nN decay mode. Widths listed in the second and third 
columns are full widths. It can be seen that we get distinct resonance signals 
even for resonances with a branching fraction as small as 10-20 % to the nN 
channel. 

Next, we move on to the analysis of the / = 3/2 partial waves in ttN scat- 
tering (Fig. 5). Many conventional pole positions of the energy-dependent 
T-matrices occur in regions of negative time delay. However, the time delay 
plots made from fits to the single energy values of the amplitude do show 
small peaks due to the A resonances with a small branching fraction to the 
ttN decay mode [37]. With the exception of P31 where we do not find any 
positive region, the overall picture is similar to the case of the 7 = 1/2 reso- 
nances. It is however fair to say that in the cases of D33, F35 and F37 partial 
waves, the positive peaks are really too tiny to be conclusive. A comparison 
of the time delay peaks in the 1 = 3/2 partial waves (evaluated using the 
FA02 solution) with the Summary Table values and pole positions of FA02 
and Speed Plots is given in Table 3. 
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Figure 5: (a) Real (dashed lines) and imaginary (solid lines) parts of the 
T-matrix solutions FA02 [56] for isospin, / = 3/2, partial waves, (b) The 
time delay evaluated using the T-matrix solutions FA02 (shown in (a)) and 
an earlier version SM95 by the same group. 



18 



Table 3. Comparison of A resonance parameters from time delay evaluated 
using the FA02 T-matrix solution, with the pole positions [62] of the same 
T-matrix, Summary Table values and Speed Plot poles (P = E - ir/2). 



L27,2J 

Mass (^rUll - 
width) 


Speed plots 

role [r ) 
Re P [-2 Im P] 


FA02 

role [r ) 
Re P [-2 Im P] 


Branching 

to vriN 
decay mode 


Time delay 

Peak [Partial 
- width] 


^33 

1232 (120) 


1 onn [1 nnl 
izuy [iUUJ 


1011 ri m 1 
izii [iUiJ 


^ on 0/ 

> yy /o 


iziU [iUoJ 


33 

1600 (350) 


iOOU [ - J 




iU - zo /o 




'J3I 

1620 (150) 


1608 [116] 


1594 [114] 


20 - 30 % 


1550 [54] 


^33 

1700 (300) 


1651 [159] 


1633 [254] 


10 - 20 % 


t 


Q*** 
'J3I 

1900 (200) 


1780 [ - ] 




10 - 30 % 




1905 (350) 


1829 [303] 


1832 [239] 


5 - 15 % 


t 


13**** 

^31 

1910 (250) 


1874 [283] 


1781 [493] 


15 - 30 % 




33 

1920 (200) 


1900 [ - ] 








'-'35 

1930 (350) 


1850 [180] 


1918 [277] 


10 - 20 % 


1800 [157] 


37 

1950 (300) 


1878 [230] 


1871 [236] 


35 - 40 % 


t 



t Resonance signal not very clear 
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4 Time delay and Breit-Wigner amplitudes 



In 2.2, we evaluated the time delay from a simple Breit-Wigner (BW) ampli- 
tude, assuming that the branching ratio for the clastic channel under consid- 
eration is 1. We discuss this topic at this point again as we wish to interpret 
the results of our time delay analysis using T-matrix solutions which involve 
BW-like functions. 

Generahzing the case of the BW amphtude by considering an elastic chan- 
nel i with a branching ratio, Br = Fj/F, smaller than one, we have 

r = ^ (11) 

Er-E- iV/2 ^ ' 

from which, after taking into account that the ^'-matrix is 5* = 1 + 2iT = 
■qe^^^ , the phase shift 5 can be calculated to be 



e- 1 -1 

6 — - tan 



r.(^R - E) 



2 L(^^-^)2 + rV4-r,F/2 

The derivative of the phase shift taken at £^ = Er is. 



(12) 



d5 
dE 



' * (13) 



E=En r {Br - 1/2) 



Prom Eq. 13, we see that if < 1/2, time delay at resonance is negative. 
The negative region around Er is a local negative minimum accompanied 

on two sides by local positive maxima. We refer to this phenomenon as 
'two-horn' structure. We can either say that the time delay concept loses 
its meaning in resonance physics in channels with a branching ratio less 
than 1/2, or, the simple Breit-Wigner is not an adequate description for a 
resonant amplitude; although for many purposes other than time delay, a 
reasonable approximation. Indeed, Eq. (11) lacks the threshold factor and 
energy dependent widths. The threshold should at least be implemented 
correctly by F, — >• Ti[q{E)/q{ER)]'^^~^^ (see e.g. [63] and references therein), 
where q is the momentum of one of the initial particles in the CM frame. This 
Fj is not unique and can be modified [54]. Moreover, as noted in [54] the rest 
of the energy dependence can take many different forms (see [64] for the latest 
collection of such resonant amplitudes). Hence, any argument in connection 
with time delay relating the latter with a Breit-Wigner should not be based 
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on (11). This is exactly what we have done by computing the time delay 
from T-matrix solutions which are based on a more sophisticated BW form 
[54] . Certainly in some cases, an almost 'two- horn' structure is visible and is 
reminiscent of a BW. However, not always should this fact be considered as 
a drawback for the time delay concept. Indeed, in spite of an almost 'two- 
horn' structure for P13, we find a mass of 1585 McV as compared to the pole 
value of 1584 MeV. The partial wave, in spite of the almost 'two-horn' 
structure, displays two bona fide resonances at their expected values. This 
shows that a more sophisticated BW can indeed account for resonances in 
the time delay. The Z^is, 6*31 and D35 resonances are in agreement with the 
values found in [54]. We do find the Pii(1710) in spite of a branching ratio 
of 10 — 20%. Hence, this again demonstrates that arguments based on (11) 
are not at all stringent. 

There are certainly problematic cases such as the D33, F37 and F35. Dis- 
missing the concept of time delay on account of these cases would amount to 
the same as dismissing the solutions of [54] which also fail to find several im- 
portant resonances quoted in PDG (in addition to the fact that the resonance 
parameters do not always agree with the mean PDG value). This is clearly 
unacceptable as it is rather the rule than exception that different analyses in 
hadronic resonances yield different results. The fact that time delay is inde- 
cisive in certain cases could be due to the BW used in the paramctrization of 
the T-matrix solutions. Although much better than (11) (as already proved 
by the time delay method itself), it might still not be the most general and 
suitable form for broad resonances [65] . In the i — > j channel, the numerator 
of the BW amplitude gets replaced by QCj with q = ^JVi/V. It was noted 
in [66] in connection with vrvr scattering, that in order to obtain the correct 
S-matrix, the q's should be complex. This should apply equally well to the 
baryon resonances (a similar argument can be found in [12]). It is worth 
noting that in our time delay analysis of meson- meson scattering [39] , in all 
cases with Sr ^ 1/2 for the elastic channel we found positive peaks (and no 
'two-horn' structures). These cases are not isolated as there are six of them. 
Hence in contrast to the simple theoretical example of an oversimplified BW, 
in reality (at least for the meson-meson case) time delay works. We therefore 
suspect that an improvement on the BW in the baryon case will also lead to 
an improvement of the corresponding time delay results. 

Finally, we also examine the influence of the non-resonant background 
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using the BW parametrization for coupled channels as in [67]. Assuming a 
diagonal background, the ^'-matrix for elastic scattering [i = j) is given as, 



1 - 2i- 



E-'-El + iERTiE) 
where r]i is the background phase. In this case. 



d5_ 
dE 



Br 



+ 



drji 



E=En r{ER){Br - 1/2) dE 
which generalizes Eq. (13). 



E=E 



(14) 



(15) 



R 



5 Summary 

Though the resonances in nN scattering belong to one of the oldest topics of 
particle and nuclear physics, their study is far from being complete. This can 
be seen alone from their classification into four-, three-, two-, and one-star 
resonances according to the status of being well or less established. Disagree- 
ments in the extractions of resonance parameters are often encountered and 
theoretical model calculations sometimes predict new resonances. 

Motivated by statements in literature that a positive time delay is a 
necessary condition to confirm a resonance, in the present work we have put 
the time delay method to test and have presented a systematic survey of time 
delay plots for the nN resonances. In case of a clear signal, the resonance 
parameters could be extracted from the time delay plots. When At{E) was 
calculated from the T-matrix solutions, we did not always find resonances 
at the energies corresponding to the poles of the T-matrix. This might be 
due to the model dependence of the T-matrix solutions or the non-resonant 
background in them. The calculation of At{E) directly from data via phase 
shifts also characterized several established resonances. The results indicate 
that in some cases a known resonance could actually be a mixture of two 
neighbouring resonances. Detailed fits to the structure in the data revealed 
new resonances. For example, in the much talked about Su partial wave 
[5, 59, 68] we found some new resonances with some of them in agreement 
with the predictions in literature [59]. We believe that given more precise 
data, the time delay approach to resonances would be a very useful tool to 
characterize resonances. 
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APPENDIX 

We briefly discuss two concepts related to time delay. These were not 
used in the present work, but we include them to avoid confusion among 
different concepts. 

Time delay is closely related to the lifetime matrix Q, defined by Smith 
[18]. This Q is related to the scattering matrix S as, Q = —ihSdSydE. 
The average time delay can be defined as a weighted average of the delay 
times Atij and is the same as Qa. Since the particle has probability \Sij\'^ of 
emerging in the j*'' channel, the average time delay for a particle injected in 
the i*'* channel is given as, 

{Ati)av = Y.S*jS^jAUj = Re 
j 

Strictly speaking, equations (1), (6) and (10) represent a time "delay" 
due to interaction. One can also derive an expression [16] (quoted here for 
J = 0): 

^/ N ./w f d8 ^ sin[2(5 + /c'a)] 1 
T{a)^%^h dE'\A{E',E)\H2— + 2a ^-^^ ^| , 

which is interpreted as the time spent within the spherical interaction region 
of radius a in the presence of interaction. Since the above equation can be 
rewritten as an integral over a probability, T(a) is always positive definite in 
contrast to (1), (6) and (10). 

We also note that yet another discussion of time delay in classical and 
quantum scattering theory can be found in [31]. 

Acknowledgements: The authors wish to thank R. A. Arndt for useful 
discussions and I. Strakovsky for providing the pole values of their T-matrix 
solutions. 
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